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Abstract. We have calculated optical spectra of hydrogen-rich (DA) white dwarfs with 
magnetic field strengths between 1 MG and 1000 MG for temperatures between 7000 K and 
50000 K. Through a least-squares minimization scheme with an evolutionary algorithm, we 
have analyzed the spectra of 114 magnetic DAs from the SDSS (95 previously published plus 
14 newly discovered within SDSS, and five discovered by SEGUE). Since we were limited to a 
single spectrum for each object we used only centered magnetic dipoles or dipoles which were 
shifted along the magnetic dipole axis. We also statistically investigated the distribution of 
magnetic-field strengths and geometries of our sample. 



1. Introduction 

The Sloan Digital Sky Survey (SDSS) project covering 10^ deg^ of the sky has discovered 
thousands of new white dwarfs, among them 105 with magnetic fields (MWDs) [1-3] .White 
dwarfs with magnetic fields between 10^ and 10^ G are thought to represent more than 10% of 
the total population of white dwarfs [4]. The SDSS Data Release 3 (DR3, ^ 4200 deg^) increased 
the number of known magnetic white dwarfs from 65 [5; 6] to 170 [7]. 

In this work we present the re-analysis of the 95 DA MWDs discovered by Schmidt et al. [2] 
and Vanlandingham et al. [3], plus the analysis of 19 additional objects from SDSS up to DR6^ 
(9583 deg2). 

Schmidt et al. [2] and Vanlandingham et al. [3] determined the field strengths and the 
inclinations of magnetic dipoles by comparing visually the observed spectra with model spectra 
calculated using a simplified radiation transfer code [8]. Their analyses resulted in dipolar field 
strengths for the SDSS MWDs between 1.5 MG and ~ 1000 MG. From their sample of 105 
MWDs, 97 were classified as DAs. 

2. SDSS data 

The spectroscopic targets of the SDSS are selected based on color criteria optimized for galaxies 
and quasars. However, follow-up spectroscopy of many stars is also performed with the twin 
dual beam spectrographs (3900-6200 and 5800-9200 A, X/dX - 1800), in particular of blue 



^ http://www.sdss.org/dr6/ 



objects like white dwarfs and hot subdwarfs [9; 10]. Since the energy distribution of strongly 
magnetic white dwarfs can differ from nonmagnetic ones, MWDs are not only found in the 
SDSS color categories for white dwarfs or blue horizontal-branch stars, but may also fall into 
the color categories for quasars (QSOs), "serendipitous blue objects'", and hot subdwarfs. Based 
on their colors, objects are assigned to fibers for spectroscopic investigation (for spectroscopic 
target selection, see [11]). 

In order to identify magnetic white dwarfs among the white dwarfs or other categories, 
different techniques were used: Prom the sample of white dwarfs, selected by color cuts, Gansicke 
ct al. [1] and Schmidt et al. [2] used visual inspection. Vanlandingham et al. [3] inspected only 
those objects visually for which bad were obtained by their autofit process. This procedure 
selects in particular white dwarfs with magnetic fields above 3MG and misses objects with 
weaker magnetic fields. 

In addition to these already known objects, we have analyzed data of nineteen additional 
SDSS objects discovered to be magnetic due to suspicious radial velocity measurements by 
the HYPERMUCHFUSS project (see, Geier et al, these proceedings and Tillich et al, these 
proceedings). The one-dimensional spectra which we used in this work were generated by SDSS's 
spectroscopic pipeline spectro2d and downloaded from Data Archive Server. 

3. Analysis 

In order to increase efficiency, we pre-computed a three-dimensional grid of Stokes / and V 
model spectra with effective temperature 7000 K < T^s < 50000 K, magnetic field strength 
1 MG < B < 1 GG, and 18 different directions ip relative to the line of sight as the independent 
variables (9 entries, equally spaced in cosV") using the radiative transfer code for magnetized 
white dwarf atmospheres [sec 13; 14]. All spectra are calculated for a surface gravity of log g = 8. 
Since no polarization information is available for the SDSS, our analysis is limited to the flux 
spectra (Stokes parameter /). 

The magnetic held geometry of the magnetic white dwarfs was determined with a modified 
version of the code developed by Euchner et al. [15]. This code calculates the total flux (and 
circular polarization) spectra for an arbitrary magnetic field topology by adding up appropriately 
weighted model spectra for a large number of surface elements. Complex magnetic field 
geometries are accounted for by a multipole expansion of the scalar magnetic potential, but 
in our paper only centered and offset dipoles were considered. The observed spectra are fitted 
using an evolutionary algorithm [16] with a least-squares quality function. 

Due to the lack of a consistent theory for Zeeman and Stark broadening, the latter is only 
taken into account by a rather simple approximation [13] so that systematic uncertainties, 
particularly in the low-field regime (< 5MG), are unavoidable. Consequently, effective 
temperatures and surface gravities derived from fitting the Balmer lines alone are less reliable 
than in the case of non-magnetic white dwarfs. This may also result in disagreements with 
temperatures estimates derived from the continuum slope. 

Our fitting procedure had to adjust three or four free parameters: the magnetic dipole field 
strength B^, the effective temperature T^s, the inclination of the dipole axis i, and an offset 
along the magnetic axis Zofltj if offset dipoles are used. Tcfj needed to be independently calculated 
for our purposes. For the 95 MWDs analyzed, we used the literature values for T^s which were 
determined by comparison to the theoretical non-magnetic DA colors in the u~ g vs g — r plane 
[2; 3]. The temperature of the new MWDs presented in Table 1 were estimated by the synthetic 
SDSS color-color diagrams by Holberg & Bergeron [17]. 

All fits resulted in reduced values between 0.8 and 3.0 except for some high-field objects. 
In Pig. 1 we show fits of 6 MWD spectra as an example. 
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Figure 1. Representative sample of fits of observed spectra of MWDs from the SDSS to centered 
magnetic dipoles with a polar field strength Bp (left) and dipoles shifted by z^fi stellar radii along 
the dipole axis (right). 
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Table 1. Best fit parameter values for centered dipole and ofi^set dipole models for the objects 
discovered by HYPERMUCHFUSS. 

^Inclination errors are estimated as large, for explanation see Sec. 4 



4. Discussion 

Three objects analyzed by Schmidt et al. [2] and Vanlandingham et al. [3] are omitted in this 
work. SDSSJ05959. 56+433521. 3, Glll-49 was listed by Schmidt et al. [2] as a magnetic DA, 
but is a Carbon-rich (DC) MWD [18]; SDSSJ084716.21+U84204 is a DAH+DB binary, hence 
we were unable to analyze these two objects with our code. Finally, neither the coordinates 
nor fiber identifier of SDSSJ23460544 +385337.7 mentioned by Vanlandingham et al. [3] was 
accessible through the SDSS database. 

Emission lines were found in SDSSJ102220.69+272539.8 and SDSSJ102239.06+194904.3 
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Figure 3. Histogram of magnetic white 
dwarfs in equal intervals of logB. Gray 
columns represent the number of all DA 
MWDs and black shades represent the 
the contribution of SDSS to DA MWDs. 
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Figure 4. Scatter plot of dipole magnetic field 
value vs. temperature in this work. The random 
distribution of field strengths with respect to the 
age indicator temperature is consistent with a long 
decay timescale of MWDs. Gray dots indicate the 
objects from HYPERMUCHFUSS. 



(the latter is shown in Fig. 1), very similar to SDSSJ121209. 31+013627.7 which could indicate 
that these objects may be EF Eri like, magnetic cataclysmic variables with a brown dwarf 
companion [19]. 

The spectra of the high-field objects SDSSJ22474141+145638.8 and SDSSJ101805.04- 
+01123.5 (PG 1015+014, shown in Fig. 1) do not fit particularly well. At higher field strengths 
(> 50 MG) the spectra become very sensitive to the details of the magnetic field geometry, as 
was demonstrated by Euchncr ct al. [15, 20, 21]. The deviations of the observed spectra from our 
theoretical spectra assuming (offset) dipole models hint therefore to a magnetic field geometry 
that is more complex than a shifted dipole. A more comprehensive analysis of PG 1015+014 [21] 
showed that individually tilted and off-centered zonal multipole components with field strengths 
between 50-90 MG are needed to represent the global magnetic field, consistent with our analysis. 
On the other side, some high-field objects in our sample, like SDSSJ135141. 13+541947 (Fig. 1) 
are well fitted. 

SDSS has proven to be extremely important in increasing the total number of MWDs. Out of 
8000 WDs (DR4; [22]), a total 139 MWDs were discovered [DR6; 23], stiU a very low percentage 
compared to the sample of MWDs in the solar neighborhood (13% ± 4%; [7]). Note, however, 
that no systematic search for MWDs in DR4-DR6 was performed yet. SDSS is designed to 
investigate galaxies and quasars, thus stars are not its primary targets. Formerly the WD 
acquisition efficiency of the SDSS was partly defined by the probability of WD spectra to be 
target labelled as QSO. With each new data release, target selection criteria has improved in 
separating the QSO from non-QSO targets. Unfortunately the improvement in extragalatic 
target selection impairs SDSSs capability to detect serendipitous MWDs. This was examplified 
by the inability to recover most of the formerly known MWDs in regions of the sky covered in 
DR2-3 by Vanlandingham et al. [3]. 

Overall, SDSS has nearly tripled the number of MWDs, conversely the completeness of the 
total MWD population is affected significantly by the SDSS biases because of this high impact of 
SDSS. The most important selection effect arises from the priority selection of the spectroscopic 
targeting, which is determined from color categories of each object [11]. Due to unconventional 
spectrum and absorption features, MWDs with high field strengths have a tendency to be 
labeled as interesting objects (e.g., BLUE_SERENDIPITY, which is an object in undocumented 
coordinates with blue colors). This contributes as a positive bias for MWDs with high magnetic 
field strengths {B > 100 MG). On the other side, seven objects from the HYPERMUCHFUSS 



contributed to the B < 3.5 MG range . This hmit was noted to be a negative bias by [3], due 
to the inscnsitivity of the autofit process. Since spectra from objects with such low fields only 
slightly differ from the non-magnetic cases, they were not found by the automatic fitting process 
and therefore no visual inspection was triggered. Since some HYPERMUCHFUSS objects were 
found in DRl to DR3, the radial velocity method has proven to be complementary in finding 
some of the missing low-field MWDs. 

These selection effects do not apply to the second phase of the SDSS, Sloan Extension for 
Galactic Understanding and Exploration (SEGUE), which is targeting the stellar population 
specifically. This includes the white dwarfs and blue horizontal-branch stars^ (that are also in 
MWD color criteria). In our sample, five out of the nineteen new MWDs are from SEGUE. 

5. Conclusion 

Overall, our results are consistent with the former analyses of MWDs (see Fig. 3), which 
shows that simple atmosphere models with preassumed dipole magnetic values are good 
approximations for these objects, if only single-phase spectroscopy without polarization 
information is considered. Our method is able to account for a diversity in low magnetic field 
strengths (B < 20 MG). This is apparent in the deviation of our best-fit magnetic field strength 
values from literature values in the horizontal ~2 MG axis of Figure 3. In many cases offset dipole 
models resulted in significantly better fits than the models with centered dipoles, however, at 
smaller field strength this can be in part due to the simplifications in the treatment of Stark 
broadening. 

The distribution of the magnetic field strengths of the MWDs from the Schmidt et al. [2] 
sample had a maximum around ~5-30 MG. We have updated the values of DA white dwarfs 
and created a histogram of all known magnetic DA WDs (Fig. 3) and added results for the new 
objects. The magnetic field strengths of non-SDSS MWDs are from the compilation by Kawka 
et al. [7]. In spite of the SDSS biases, the same ~5-30 MG peak as in Schmidt et al. [2] is 
apparent in Fig. 3. 

The dipole magnetic field ohmic decay timescale is lO^'^ yr. Even the higher multipoles can 
live for such a long period of time (see, e.g. [24]). Therefore, no significant correlation between 
temperature and magnetic field strength is expected if temperature is assumed as an indicator 
of age (Fig. 4) . This lack of correlation supports the fossil ancestry of these fields inherited from 
earlier stages of stellar evolution. 

High field MWDs are thought to be remnants of magnetic Ap and Bp stars. If fiux 
conservation is assumed, the distribution of the polar field strengths of high field MWDs should 
be largest in the interval 50-500 MG. In our sample, objects with magnetic field strengths 
lower than 50 MG are more numerous than the objects with higher magnetic field strengths 
(see Fig. 3), part of this effect is due to our biases (see Section 2). Nevertheless it is consistent 
with previous results and supports the hypothesis that magnetic fossil fields from Ap/Bp stars 
only are not sufficient to produce high field MWDs [25]. Auricre et al. [26] argued that dipole 
magnetic field strengths of magnetic Ap/Bp stars have a "magnetic threshold" due to large scale 
stability conditions, and this results in a steep decrease in the number of magnetic Ap/Bp stars 
below polar magnetic fields of 300 G. 

A possible progenitor population for MWDs with dipolar field strengths below 50 MG is the 
currently unobserved population of A and B stars with magnetic field strengths of 10- 100 G. 
Wickramasinghc & Fcrrario [25] suggested that if 40% of A/B stars have magnetism, this 
would be sufficient to explain the observed distribution of MWDs. However, the existence of 
this population seems to be highly unlikely since investigations of Shorlin et al. [27] and Bagnulo 
et al. [28] for magnetism in this population yielded null results, for median errors of 15 - 50 G 

^ http://segue.uchicago.edu/targetsel.html#v4.2 



and 80 G respectively. Another candidate group is the yet undetected magnetic F stars for these 
MWDs with lower field strengths [2]. But this conclusion is strongly affected by SDSS MWD 
discovery biases. 
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